Constraints on charged Higgs bosons from 
Df ^jL^u and Df r^u 



A.G. Akeroyd^l and F. MahmoudiH 



a; Department of Physics, National Central University, Jhongli, Taiwan 320 
b: Laboratoire de Physique Corpusculaire de Clermont-Ferrand (LPC), Universite 
Blaise Pascal, CNRS/IN2P3, 63177 Aubiere Cedex, France 



Abstract 

The decays — > /-i^f^ and r^Vr have traditionahy been used to mea- 

sure the Df meson decay constant //j^. Recent measurements at CLEO-c and the 
B factories suggest a branching ratio for both decays somewhat higher than the 
Standard Model prediction using f^^ from unquenched lattice calculations. The 
charged Higgs boson {H^) in the Two Higgs Doublet Model (Type H) would also 
mediate these decays, but any sizeable contribution from can only suppress 
the branching ratios and consequently is now slightly disfavoured. It is shown that 
constraints on the parameters tan /? and mij± from such decays can be competitive 
with and complementary to analogous constraints derived from the leptonic meson 
decays t^Vj- and Ai^t'^, especially if lattice calculations eventually 

prefer fo, < 250 MeV. 



1 Introduction 



In the Standard Model (SM) the purely leptonic decays of charged pseudoscalar mesons 
(P^ — ^ i^u) proceed via annihilation in the s-channel to a. W boson. The decay rates 
are proportional to the lepton mass which arises from the chirality flip of the lepton 
required to conserve angular momentum. Such decays have traditionally been used to 
measure the decay constants of the pseudoscalar mesons, and thus provide an important 
test of lattice Quantum Chromodynamics (QCD) calculations. The unprecedented data 
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samples provided by the B factories have enabled the decay r^z/ to be observed 

P, El El H] despite its relatively small branching ratio (BR) and challenging signature. 
The decays ^ ^i^v Df ^ and Df r±i/ [3 El [IHl [II] have 

been measured with ~ 10% precision at CLEO-c, and the recently commenced BES-III 
experiment will provide improved measurements [T^fTB] . Moreover, the decay Df ^ fi^u 
has been measured at the B factories [HI [I5] with a precision slightly less than that of 
CLEO-c. 

These purely leptonic decays are sensitive to charged Higgs boson (H^) at the tree level 
and thus provide valuable probes of such particles which are complementary to constraints 
provided by loop-induced decays (e.g., b — > 57). Importantly, in supersymmetric (SUSY) 
models the loop-induced decays are particularly sensitive to the sparticle spectrum and the 
assumptions made for the SUSY breaking sector, and thus the purely leptonic decays offer 
a more model-independent probe of parameters in the Higgs sector. The measurement of 
B^ T^u [HE], although in rough agreement with the SM prediction, provides important 
constraints on the mass and couplings of from the Two Higgs Doublet Model (Type II) 
[T6] . with some dependence on SUSY parameters [171 [TBI [19l [20l EH ESj which enters at 
higher orders in perturbation theory. Improved measurements of B^ — ^ t^u (and first 
observation of B^ — > n^i') are thus certainly desirable and such decays play a prominent 
role in the physics case for a high luminosity flavour factory [231 EH ESI ESI EZl EHl EHl ESI 
EI]. 

The decay rates of fi^u and Df fi^u/T^u have traditionally been considered 

as a robust test of lattice QCD calculations of the decay constants Jd and Jd^ because 
significant New Physics effects were thought to be unlikely. However, the effect of 
on Df — > i^^u /t^u was pointed out to be non-negligible in [ISl E21 ES] and a numerical 
study [Ml ES] showed that the magnitude of the contribution can be comparable 
to both the CLEO-c error and the lattice QCD error for the decay constants. SUSY 
particles from models with i?-parity violating interactions can also contribute to both 
D^ — > /i^i/ and Df fi^u/r^u [3S]. Recently, there has been growing interest in the 
effect of New Physics particles on Df /i^z//r^z/ [37] because the current world average 
of the experimental measurements of their branching ratios is somewhat higher than the 
SM rate using unquenched lattice QCD calculations of fn^ [SHIES]- Such a discrepancy 
cannot be explained by of the popular Two Higgs Doublet Model (Type II) since any 
sizeable contribution (which arises from a relatively large strange quark Yukawa coupling) 
can only suppress the decay rate [ISl E2l EH ES]- In this paper we consider in detail the 
effect of on Df —>■ fi^u/T^u using the program Superlso [ISl HI]. It is shown that 
constraints derived on the plane [tan/3, mj|^±] can be competitive with the constraints 
obtained from the analogous leptonic decays B"^ r^v and n^u. In particular, 

we compare the potential of Df — > h^u/t^u and — > jji^p [121 ESI HI] to probe the 
plane [tan/3, mj:^±]. 

The paper is organized as follows: in section 2 the theoretical formalism for the decays 
Df ^i^v/t^v and — > jji^u is introduced; numerical results are contained in section 
3 with conclusions in section 4. 
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2 Effect of H on leptonic decays of Df and K 

In this section we review the theoretical formahsm for the effect of on Df — > /i^z/, 
Df —* t'^u and — > n^u, and summarize the experimental situation. 



2.1 The decays Df ^ /i^u and Df ^ t^u 

In the SM the purely leptonic decays Df i^U£ proceed via annihilation of the heavy 
meson into W*. Singly charged Higgs bosons, which arise in any extension of the SM 
with at least two SU{2)l x U{1)y Higgs doublets with hypercharge Y = 1, would also 
contribute to these decays [I6]. The tree-level partial width is given by (where i = e, n 
or r): 

T{Df tu,) = ^flm]Mr,^ [l - ^) Vc/-. , (1) 

where in the Two Higgs Doublet Model (2HDM) of Type II [ISl El ESI EZl EH] : 

MdA^ f m,tan2/3 
rric 



nic + msj \mH+ J \ l + eotan/5 



(2) 



Here mc and are the masses of the charm and strange quarks respectively, is 
the mass of the charged Higgs boson, M/)^ is the mass of the Df meson, tan/? = ^2/^1 
where Vi and V2 are the vacuum expectation values of the two scalar doublets, Vcs is a 
Cabibbo-Kobayashi-Maskawa (CKM) matrix element, mi is the lepton mass and Gi? is 
the Fermi constant. In the 2HDM (Type II) each fermion receives mass from one vacuum 
expectation value {yi or V2) at tree level and consequently eo = 0. In the 2HDM (Type 
HI) each fermion receives a mass from both vi and V2 and the term eo 7^ must be 
included. The Yukawa couplings of the Minimal Supersymmetric SM (MSSM) take the 
form of a 2HDM (Type II) at tree level, but at higher orders the structure becomes of the 
type 2HDM (Type HI) in which eo is a function of SUSY parameters [15], HEl HTJ HE] and 
|eo| can reach values of order 0.01. 

The above formula Eq. ([2]) was derived in [16] for eo = and neglecting the term 
in {rric — nisia.v? (3). This term originates from the charm quark Yukawa coupling 
and increases r<, (and was included in [32]) but its magnitude is numerically tiny and can 
be neglected given the current experimental and theoretical errors. The term m<j tan^ /3, 
which originates from the strange quark Yukawa coupling, can give rise to a non-negligible 
suppression of for large values of tan (5. Both [16] and [32] neglected when it appears 
as (mc + m^) and this correction was first included in [2S]- The above complete expression 
for the case of eo = is taken from [STJ EH]- The eo correction can be non-negligible for 
large tan j3 and should be included in a dedicated study of Ts in the context of the MSSM. 
The parameterization of the eo correction in Eq. ([21) is valid in MSSM models with minimal 
flavour violation. We neglect to write explicitly the analogous eo correction for the rric 
term since this correction does not have the tan (3 enhancement factor, and the term 
is negligible anyway. A first quantitative study of the magnitude of (setting eo = 0) 
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was performed in [Ml [35] for the tan^ /3 term only and it was shown that could be 
suppressed by an amount comparable to both the CLEO-c experimental precision and 
the theoretical error in the lattice calculations of fn^. Note that the magnitude of the 
H"^ contribution depends on the ratio of quark masses ms/{mc + ttLs), and an analogous 
uncertainty is not present for the contribution to the decay — > i^u. 

There are various unquenched lattice calculations of Jd^ and the current situation is 
summarized in [38] and updated in [9], HH]. The value with the smallest quoted error is 
Jd^ = 241 ± 3 MeV [50] which employs staggered fermions. Another calculation with 
staggered fermions gives f^^ = 249 ±11 MeV [ill EI]. Other calculations use different 
lattice techniques and have larger errors than that of [HO], although usually with a central 
value less than 250 MeV (e.g., 248±3±8 MeV [52]). On the experimental side, Df ^ r^i^ 
has been measured at CLEO-c for two decay modes of r [9], [H] and Df has 
been measured at CLEO-c [9], BELLE [H] and BABAR [15]. Taking an average of 
these measurements (excluding the BABAR result which has an additional error from 
normalizing to the branching ratio of Df 0°7r^) results in an unexpectedly high value 
of = 261 ±7 MeV (derived in [9] and is dominated by the average of the three CLEO-c 
measurements in [9] and [H]). A previous average gave = 273 ± 10 MeV [381 139] . 
The sizeable difference between this world average measurement of fn^ and the lattice 
QCD value with smallest error {fos = 241 ± 3 MeV [SU]) has brought attention to the 
effect of New Physics particles on these decays [H^. As discussed above, if^ in the 2HDM 
(Model II) can only give a sizeable suppression of the branching ratio for Df yU^i^, t^i^, 
and not a sizeable enhancement. Candidate models which can enhance the BRs sufficiently 
are -R-parity violating models [Ml EH |53] or models with leptoquarks [371 El]- We note 
that the scenario of a charged Higgs boson from an unconventional 2HDM with a large 
charm quark Yukawa coupling [27] (which is never possible in standard 2HDMs in which 
the top and charm quark receive mass from the same vacuum expectation value) was also 
available to accommodate the data for BR(D^ — > h^u/t^u). However, such a would 
enhance BK{D^ — > fi^i') by the same amount [55] and now appears to be disfavoured by 
the final CLEO-c measurement of BR(D^ — s> ^^u) and fu [6], which is in good agreement 
with the lattice calculation of fn in [30] . 

In [nS] two lattice QCD calculations [HOI EI] of fn^ were combined with the value 
obtained from QCD spectral sum rules in order to give an average value of Jd^- Moreover, 
an analogous average of calculations of fo and fos/ fn were multiplied together to give 
an independent evaluation of fn^- The above procedure gives the value fr>^ = 240 ± 7 
MeV (i.e. very similar to fo^ = 241 ± 3 MeV of ^U])- The consequences of relaxing the 
assumption of unitarity of the CKM matrix when extracting the measurement of from 
BK{Df — > fi^u/T^u) were also discussed. In the models we study in this work (2HDM 
and MSSM) the CKM matrix is unitary and \Vcs\ ~ \Vud\ with very small error. As 
discussed in [56], this scenario leads to the largest discrepancy between the experimental 
and theoretical values of fu^. 

In order to make conclusive statements about the possibility of New Physics enhancing 
Df fi^u/r^u it is evident that a detailed discussion of the errors for the value of fo^ 
evaluated in [50] is required, as well as a consensus in the lattice community concerning 
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the magnitude of the error. Our view is that the current situation of the world average 
measurement of fo^ being higher than most lattice calculations of Jd^ is provocative and 
certainly merits attention because more precise measurements of Df — > fi^u/T^u will be 
possible at the recently commenced BES-lII experiment [12]. In addition, high luminosity 
flavour factories are also being discussed [23] which have the option of running at the charm 
threshold [57] with data samples even larger than those anticipated at BES-III. In this 
work our aim is not to explain the currently higher than expected BIi{Df yU^z//r^z/) 
by invoking speciflc New Physics models. Indeed, one cannot predict the preferred values 
for the lattice calculations of fn^ although one expects a gradual decrease of the error. 
The current discrepancy between the experimental and lattice QCD values for f£,^ might 
persist, worsen or disappear entirely. Instead, we wish to emphasize that BR(D^ — ^ 
fi^u/r^i') should be considered seriously as a process which potentially offers constraints 
on models containing a with the popular structure of the 2HDM (Type II), which 
includes the MSSM. 

In [3HI El] the SM prediction for BR(D^ — + ii^u^/t^Ut^) using the lattice calculation 
with smallest error {fos = 241 ± 3 MeV [50]) was compared with the experimental world 
average at the time for Jd^ (273 ± 10 MeV). It was concluded that any contribution 
is disfavoured at more than 3a. Such a result shows the potentially important role of the 
decays Df — > fi^u/r^u as a probe of H^. However, it is also clear that constraints on 
derived from these decays are strongly dependent on the central value and error for 
the lattice calculation of fo^ and so should be interpreted with caution at present. In this 
work we aim to quantify the constraints on the plane [tan (3, mH±] from Df — *• ^"^v^J t^Ut 
with various assumptions for the central value and error of the lattice QCD calculation of 
ff)^ , and compare the constraints to those obtained from other decays which are sensitive 
to 
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2.2 The decay K 

A decay which has many similarities with Df fi^u and Df t^v is — > /i^i/. This 
decay can also be mediated by of the 2HDM (Type II) at tree level [TB] with a scale 
factor similar to that in Eq. ([2]). In order to reduce the theoretical uncertainties from the 
decay constant /i^, the ratio of partial widths is usually considered: 



r(K± 



v., 



ud 



fK^K f 1 - mj/mj^ 



X 



m 



K+ 



m 



tan^ (3 



1 + eo tan P 



(3) 



Here 6em = 0.0070 ± 0.0035 is a long distance electromagnetic correction factor. As 
suggested in [12], we study instead the quantity 



123 



Vus{Kn) Kd(0+ ^ 0+) 



X 



i3) 



(4) 
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Here VusiKu) refers to Vus as measured in leptonic decay of K"^ with i particles in the final 
state (two leptons and a number of pions), and similarly for Vud- Vud(0^ —>■ 0+) denotes 
Vud measured from nuclear beta decay. In the SM i?£23 = 1? while the contribution from 
in the MSSM attains the simple form (which can be compared to Eq. ([2])): 



R. 



■m 



mj^+ \ nis J 1 + eo tan /3 



(5) 



Using rrid/ms = 1/20 [39], the MSSM prediction can be directly compared to the experi- 
mental value [l2] 

= 1.004 ±0.007. (6) 

In the extraction of this value, the ratio fx/fn has been fixed to the value fx/fn = 
1.189 ± 0.007 obtained from lattice QCD using staggered quarks [50]. It should be noted 
that the uncertainty thus obtained for i?^23 is most probably overly optimistic. Indeed, 
many approaches exist to determine fx/fn, and some reservations remain about staggered 
fermions despite their impressive success in agreeing with the experimental values of 
various "gold-plated" observables in flavour physics (e.g., contrasting opinions concerning 
their theoretical viability are expressed in [58] and [59] )• If, for example, the value fx/f-K = 
1.205 ±0.018 (obtained using the domain wall formulation [60]) is used instead, then i?£23 
provides no constraints on the plane [tan/3,m/^±] [H]. A recent analysis in [HI] used 
fx/ f-K = 1-19 ± 0.015 and also found very weak constraints. We therefore stress that 
the constraints obtained from Eq. ([U]), although certainly valuable due to their sensitivity 
to at tree level, should serve only as an indication of the potential of K"^ decays to 
probe the plane [tan/3, mj|/±]. It is our aim to compare such constraints in this optimum 
scenario for _R£23 (i-e. using the lattice calculations of [50]) with those derived from 
Df jM^u/r^u. We note that R121, also requires lattice QCD input for the semi-leptonic 
kaon decay form factor /_|_, while the decays Df — > ^'^v/t'^v rely on lattice QCD only 
for f£)^. The direct measurement of Vus (relevant for Rm) is much more precise than 
that for Vcs (relevant for Df ^'^v/t^v). However, Vcs can be taken as a well-measured 
parameter by assuming a unitary CKM matrix, which is the case for the models we study. 



3 Numerical Analysis 

The branching ratios for Df r^Ur and Df f^'^^i^/i have been implemented in Superlso 
[101111] following Eqs. ([I]) and ([2]). The eg correction for the second quark generation has 
also been added to take into account the SUSY effect. This effect has been neglected in 
all previous studies of BR(_D^ t^J-'t) and BR(D^ — > jJi^v^). 

To illustrate the constraining power of these observables, we consider the MSSM with 
minimal flavour violation (MFV), and in particular the Non-Universal Higgs Mass model 
(NUHM). The NUHM assumes SUSY breaking mediated by gravity and is characterized 
by a set of universal parameters at the GUT scale {mo, mi/2, ^o, tan/3, /x, M^}. The 
first four parameters in this set are the same as in CMSSM (or mSUGRA), but the GUT 
scale mass parameter universality is relaxed for the Higgs sector leading to two additional 
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parameters, fj, and Ma, which imphes in particular that the charged Higgs mass can be 
considered as a free parameter. This additional freedom makes this model attractive for 
studying the constraints on the Higgs sector. 

To investigate the NUHM parameter space we generate 50,000 random points scanning 
over the ranges tuq e [50, 2000] GeV, mi/2 G [50, 2000] GeV, Ao E [-2000, 2000] GeV, 
/i e [-2000, 2000] GeV, tua G [5, 600] GeV and tan /3 e [1, 60]. For each point we calculate 
the spectrum of SUSY particle masses and couplings using SOFTSUSY 2.0.18 [62] and we 
compute the branching fractions BR(D^ — >• t'^i't), BR(D^ — > [i^p^) and BR(K='= — >• iJi^v^) 
using Superlso v2.4 [H]. The obtained values are then compared to the experimentally 
allowed intervals while taking uncertainties in consideration. 

The PDG08 [39] combined experimental results for the branching fractions of Df 
T^u and Df — fi^u are: 

BR(Df T^u) = (6.6 ± 0.6) x lO'^ , (7) 

BR{Df fi^u) = (6.2 ± 0.6) x 10^^ . (8) 

The final CLEO-c results using 600 pb~^ have been recently released [HI [TT] and give 
smaller central values and errors: 

BR(Df ^ fx^i^) = (5.65 ± 0.45 ± 0.17) x 10"^ , (9) 

BR(D,^ T^u) = (5.62 ± 0.41 ± 0.16) x 10"^ . (10) 

For BR(D^ — > T^u) the average in Eq. ([7]) is from the measurements of [3, ESI EH ES] • 
order to have an updated world average we replace the previous CLEO-c measurement of 
[7| by the new measurement [HI [TT] to obtain: 

BR{I)f T^u) = (5.7 ± 0.4) X 10"^ . (11) 

For BR(Z)^ — i> fi^i') the constraints on the plane [tan/3,m//±] turn out to be weaker 
than those from BR(Z}^ t^u). Therefore for illustration and comparison we only 
average over the CLEO-c [H] and BELLE [Tl] measurements of BR(Z}^ f^^i^) (as 
suggested in [381 139]) which are the only measurements without an additional error from 
normalizing to the branching ratio of Df —>■ (p^ir^. Our average is: 

BR(D,± ^ /i^z/) = (5.8 ± 0.4) X 10-3 . (12) 

In our numerical analysis we employ the lattice calculation of the Df decay constant with 
smallest error [50] in order to give the SM prediction for the branching ratios: 

fo, = 241 ± 3 MeV , (13) 

BR(Df T^u) = (4.82 ± 0.14) x 10"^ , (14) 

BR(Df /i^i/) = (4.98 ± 0.15) x 10"^ . (15) 

We will comment qualitatively on the case of using a different value with a larger error. 
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We use the most up-to-date PDG08 values for the quark masses |39j : 

rus = 104:1:11 MeV rric = 1.27;[J:?I GeV . (16) 

To take into account this uncertainty we vary the ratio rris/mc in the interval [0.04,0.12] 
with ms/mc = 0.08 being the central value. 

An estimation of the total theoretical error in the decay rates (Eq. [T]) yields about 3% 
relative uncertainty. The dominant theoretical error is from the lattice evaluation of fo^, 
while the error from Vcs is very small if CKM unitarity is assumed (which is the case in 
the models which we study) . The errors from SUSY and other parameters in Eq. ([T]) yield 
less than 1% relative uncertainty. In order to undertake a slightly conservative approach 
we consider 4% relative theoretical uncertainty in our numerical analysis. Eqs. (fTTj) and 
(IT^ lead to the following allowed intervals at 95% C.L. for BR(Df t^v): 

4.8 X 10~2 < BR(Df ^ T^u) < 6.6 x 10"^ , (17) 
and for BR{T)f fi^u) 

4.9 X 10^^ < BR(Df /i^i/) < 6.7 x 10"^ , (18) 

in which both experimental and theoretical errors are included, as in the analysis of [H]. 

Fig.[T]illustrates the obtained constraints by a projection of the six dimensional NUHM 
parameter space onto the plane [tan f3,mH±]- The allowed points are displayed in the 
foreground in green while the points excluded by BR(D^ -h> t'^u) at 95% C.L are displayed 
in the background in red. The results for three choices of ms/mc are presented. A large 
part of the parameter plane is excluded even for the pessimistic case of m^/m^ = 0.04. 
The region where red and green points overlap is caused by the theoretical uncertainty 
(mainly arising from /^J and from the effect of the eo term which may take either sign 
and depends on all six NUHM parameters. 

In Fig. [2] the points excluded by BR(D^ — s> /i^i^) are shown in black for the same 
three choices of mg/m^. The obtained constraints are less severe in this case since the 
predicted SM value is closer to the central experimental value for BR(D^ f^^'^) than 
for BR(D^ T^u). 

In order to compare our results with the constraints from K"^ — > n^u, we present in 
Fig. [3] the points excluded at 95% C.L. by the Ri23 observable using the range given by 
Eq. ([H]) with no additional error, which is probably overly optimistic as discussed in [Hj 
and in section 2.2. It is clear that BR(D^ r^v) is more powerful at constraining the 
plane [tan/3,m/^±] than i?£23 even for the pessimistic case of ms/mc = 0.04 (likewise for 
BR(D^ — >■ fi^iy) if ms/mc > 0.04). Although i?£23 has less uncertainty from quark mass 
parameters than BR(D^ — ^'^v/t^v) (i.e., the theoretical uncertainty from ma/ms is 
less than that from mg/mc) we stress that the pessimistic case of mg/mc = 0.04 already 
provides competitive constraints on the plane [tan f3,mH±]- The lattice techniques dis- 
cussed in [nni promise precise calculations of and mg/mc, which would further sharpen 
potential constraints from BR(D^ ^'^v/r^v). We note that at a higher confidence level 
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°100 200 300 400 500 600 
m^, (GeV) 



Figure 1: Constraints on the plane [tan l3,mH±] from the decay Df — > t'^u for nis/mc = 
0.04, 0.08 and 0.12 from top to bottom. Red points are excluded at 95% C.L. 



9 
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m^. (GeV) 




m^, (GeV) 

Figure 2: Constraints on the plane [tan/5,m//±] from the decay Df —>■ fi^u for 
0.04, 0.08 and 0.12 from top to bottom. Black points are excluded at 95% C.L. 
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m^, (GeV) 



Figure 3: Constraints on the plane [tan/3, mj|/±] from the observable i?£23 (which involves 
— > /i^iz/vr^ — i> /i^z/ decays). Blue points are excluded at 95% CL. 



(e.g., 99.7% CL.) the constraints on the plane [tan/3, mj:^±] from BR(D^ ^^v/r^i') 
are weaker than those from i?^23 because the latter is measured more precisely. 

As pointed out in [12] , one of the virtues of -R^23 is the fact that it can probe a region 
of plane [tan/?,mij±] which cannot be probed by the decay — ^ r^v (corresponding 
to the case of the if^ induced amplitude for B"^ — > r^i/ being similar in size to that of 
the SM contribution but with opposite sign). We stress that BR(D^ t^p) can exclude 
this parameter space with greater significance than Rm-, using the lattice calculations of 
[50] for the SM prediction. Therefore both Df ^ r^z/ and Df fi^u offer constraints 
competitive with those of R£23 and complementary to those of B^ r^u. Moreover, 
we note that the experimental prospects are more favourable for BR(D^ —>■ ^'^v) and 
BR(D^ — > T^v) than for -R^23- The recently commenced BES-III experiment [12] aims to 
reduce the error of both BR(D^ f^''^^^) and BR(D^ t^v) to the level of a few percent 
(i.e. two to three times smaller than the current error) while a similar improvement for 
the precision of i?£23 seems unlikely in the same time scale. 

All the above numerical analysis has been performed with //^^ = 241 ± 3 MeV, and we 
now comment on the scenario of using a higher value of fo^ with larger errors. Repeating 
the numerical analysis with //j^ = 248 ± 3 ± 8 MeV [S2| results in almost no constraints 
on the plane [tan/5, m//±]. This conclusion also applies to constraints derived from i?£23 
using calculations of Jk/ f-n with a larger error (see section 2.2). Hence it is clear that 
strict constraints on the plane [tan/5,m/^±] from BR(D^ t^i^t/ f^'^i^fi) and -R£23 are 
strongly dependent on input from lattice QCD. However, their tree- level sensitivity to 
H"^ makes them potentially valuable probes of the MSSM and 2HDM (Type II), despite 
the above reservations. Both BR(D^ t^i^t) and BR(D^ — ^ f^^'^ti) could start to 
play a significant role in constraining the plane [tan/3, mj|/±] if future lattice calculations 
of fr)^ favour f^^ < 250 MeV. Given the ongoing interest in calculating f^^ in lattice 
QCD (summarized in [IH]) and the timely commencement of the BES-III experiment, we 
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NUHM m„=500GeV m,,2=500 GeV A„=0 u=500GeV 
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Figure 4: Regions of the plane [tan /3, m^] in the NUHM model excluded by different 
flavour observables. The constraints are superimposed in the order given in the legend. 

emphasize that these decays should not be overlooked in phenomenological studies of if^ 
in the MSSM and 2HDM (Type II). 

Finally we discuss the constraints on the plane [tan/3, mjy±] from B physics observ- 
ables. As discussed in the introduction, one-loop induced decays such as 6 — > S7 [67] 
and Bg —>■ [681 EH] can give powerful constraints on the plane [tan/5,mj:/±] (e.g. 

[211 im EOl EH [Z2]) although with a stronger dependence on SUSY parameters than the 
purely leptonic decays discussed up to now. In Fig. HI we present an example of the con- 
straints obtained by different flavour observables (both tree-level and one-loop induced) 
for fixed values of mo = mi/2 = /i = 500 GeV and Aq = 0. We show the results in 
the plane [tan/3,myi], where niA is the mass of the pseudoscalar which is very close in 
magnitude to m//± for > 200 GeV. The regions excluded by 6 — S7 observables are 
displayed in red for the isospin asymmetry and in blue for the branching ratio [721 E3] • 
The region excluded by BR(D^ — > t^u) is depicted in yellow. The green area repre- 
sents the region excluded by BR(B^ — > r^i^), the violet region by BR(Bg — > /i+/i~), the 
light blue region by Z^"*^^, and the orange area by BR(B —>■ Drz/) [H], the lat- 

ter also being sensitive at tree level to and can provide constraints competitive with 
BR(B^ T^u) [TKl [T^ fIU{ [7^ [75] . To obtain the constraints presented in this figure the 
input values of [H] are used. It is important to remember that the constraints can be 
subject to uncertainties, in particular from decay constants and CKM matrix elements. 
To obtain the constraint from BR(D^ — > t^i>) the central value ms/mc = 0.08 is used. 
Finally, the black region in the figure represents the region excluded by the direct searches 
at colliders [SH]. This figure shows that BR(D^ —* t^u) can be very competitive with the 
other decays in the context of the NUHM. 
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4 Conclusions 



Constraints on the parameter space of [tan/3, m//±] have been derived from the effect of 
singly charged Higgs bosons on the decays Df fJ'^'^^i and Df — > r^z/^ in the context 
of the Minimal Supersymmetric Standard Model. It was shown that such constraints 
can be competitive with and complementary to those derived from other flavour physics 
observables, especially if future lattice QCD calculations favour f^^ < 250 MeV (which is 
already suggested by the result of [50] )• We emphasize that such decays will be measured 
more precisely at the charm facility BES-III and hence could play an important role in 
constraining the plane [tan/3, m/^±] in the context of supersymmetric models in the future. 
The numerical analysis was performed by the publically available code Superlso [iQl SI] 
and we encourage the inclusion of Df |J'^J-'^l and Df — > r^Ur in other studies of flavour 
physics constraints on supersymmetric models. 
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